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Abstract. In this review, we provide an account of the recent progress in understanding 
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nature of disorder in these systems and its relation to transport properties. While the focus is 
primarily on theoretical and conceptual aspects, connections to experiments are also included. 
Issues such as short versus long-range disorder, localization (strong and weak), the cairier 
density dependence of the conductivity, and conductance fluctuations are considered and some 
open problems are pointed out. 
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1. Introduction 

In just the few years since its first synthesis lU E] O, graphene has become the center 
of attention of a substantial body of scientists in the fields of physics, materials science, 
chemistry, and electrical engineering. It is a fascinating material, made of a single atomic layer 
of carbon, with unique electrical, thermal, and mechanical properties. As such, it is likely to 
have a large number of practical applications. Even though the main features of its electronic 
band structure were described more than sixty years ago [4 |, only in recent years we started 
to appreciate (and understand) the wealth of phenomena in which charge carriers in graphene 
can participate. There are already a couple comprehensive reviews in the literature covering 
the fundamentals of graphene 15] |6l |2l [8] and the reader is strongly encouraged to consult 
them for acquiring a background on the subject. Here we focus instead on how disorder in 
its various forms affects the electronic transport properties of graphene at low temperatures. 
This topic has become a vast subfield in its own right. Thus, rather than provide an exhaustive 
review, in the spirit of a topical review we try to cover just enough material to bring the reader 
up-to-date with key concepts and crucial results. As a consequence, to narrow the scope, we 
leave out topics such as ac transport, quantum Hall effect, and spin transport, and the effect 
of electron-electron interactions. Due to the limited space, we also do not discuss electronic 
transport in graphene nanoribbons. 

The paper is organized as follows. In section|2]we describe the role of symmetries and 
symmetry breaking by disorder in graphene. The nature and differences among the various 
sources of disorder in graphene are discussed in section [3] In section |4] we look at transport 
properties at the neutrality or Dirac point, while in section|5]we discuss weak localization and 
conductance fluctuations caused by quantum interference. The case of doped (non neutral) 
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graphene is considered in section|6] Finally, in section|2]we draw some conclusions and point 
to a number of open questions. Although our emphasis is on theoretical aspects, connections 
to experimental results are made throughout the paper 



2. Pristine graphene and the role of symmetry 

To understand how disorder affects electronic transport properties in graphene, it is important 
to consider the ideal case of a pristine crystal Ig). Single-layer graphene is formed by 
carbon atoms disposed in a two-dimensional non-Bravais honeycomb lattice (namely, two 
intersecting hexagonal sublattices - see figure[T]i. The distance between nearest-neighbor sites 
is approximately 1.42 A, while the lattice constant is a = 2.46 A. The geometry and flatness 
of the lattice prohibit any overlap between the p- orbital of a given atom and the s, p,, and 
Pj orbitals of its neighbors. The orbitals s, p^^, and py hybridize to create the sp^ bonds that 
hold the atoms together, as well as a high-energy a band. The so-called n band is created by 
the overlap of p- orbitals and can be treated independently from other bands. In graphene, the 
n band is responsible for most of the electron conduction. Thus, at low energies or doping, 
a single-band tight-binding model including only nearest-neighbor hoping between adjacent 
sites on distinct hexagonal sublattices (which are named A and B hereafter) provides a good 
approximation for study of electronic properties in graphene. The hopping amplitude in this 
case isf «2.7eV|l5]. 
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Figure 1: Honeycomb lattice and the corresponding Brillouin zone. Lattice vectors are 
denoted by ai and 02. Empty and solid circles represent different sublattices. 

When switching to a momentum representation, in addition to the spin degree of 
freedom, one has to introduce an isospin structure in the electron wavefunction to account 
for the two sublattices. The resulting Hamiltonian can be diagonalized to yield the energy 
band El 

£±(fc) = ±fVl+4/(fc), (1) 

where 

f{k) = cos^(M/2) +cos(V3M/2)cos(/tya/2). (2) 

Neutral graphene has its Fermi energy at = 0, cutting the n band exactly at the six corners 
of the hexagonal Brillouin zone. As shown in figure [1] only two sets of nonequivalent corner 
points exist and are denoted by K and K', with K and K' points alternating along the hexagon. 
Near any of these points, say K, the two (positive and negative) branches of the n band in ([TJ 
can be approximated as 

E±{k)^±hvF\k-Kl (3) 
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where vp — ^/3at /2Ti w 10^ m/s is the Fermi velocity. This is a remarkable result: near 
E — Q (the so-called Dirac point), carriers in graphene have a linear dispersion relation, much 
like massless relativistic particles. In addition, we can introduce another isospin index to the 
wavefunction in order to differentiate contributions from the K and K' points (also known as 
the two valleys), yielding the four-dimensional spinor 



V <' / 



(4) 



(real spin indices are omitted). Close to the Dirac point, the effective Hamiltonian acting on 
this four-dimensional spinor reads 



H = vp {p ■ a) T, 
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(5) 



where the linear momentum p — h{kx,ky) is defined with respect to the corner points of the 
Brillouin zone and the Pauli matrices a = {Gx, CJy) and act on the AB and KK' subspaces, 
respectively. The eigenstates of equation (|5]l. 
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(6) 



provide a basis for solving transport problems in the presence of disorder Here, Qi^ = 
arciaaiky / kx) and the +/— signs correspond to the conduction/valence bands, respectively. 

One can use the eigenstates in equations Q to calculate the conductance in pristine 
graphene. The results are strongly dependent on geometry and edge orientation IIIOI [TTl . 
For later comparison, let us briefly comment on these results. The standard ballistic model 
considers a graphene ribbon of width W and length L connected to heavily doped graphene 
leads. In this geometry, the transverse momentum is quantized and free propagating modes 
exist at the contacts. This allows one to calculate the 5-matrix and obtain the conductance G 
using the Landauer formula. The conductivity, defined as <t = [L/W)G, has a minimum at 
the Dirac point described by 



^0- 



1 



(7) 



,1 coah^iTtnL/W) 

where go = 1/2 for metallic ribbons, and go = for semiconducting ones ifTTl . For W « L, 
<7 is dominated by go- In general, a, as given by equation (|7]l, depends on the aspect ratio 
W /L. Hence, calling it conductivity is somewhat misleading, since the conductivity is usually 
expected to be a sample specific quantity and independent of geometry. Yet, for W/L 3> 1, <7 
is dominated by evanescent modes and one finds 

This value is often (mis)quoted as the universal ballistic conductivity minimum, even though 
its derivation is only applicable to samples with a narrow aspect-ratio. The Landauer formula 
can also be used to calculate a as a function of the carrier concentration n since the latter is 
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related to the Fermi energy; n = Ep\Ep\/7t{hvp)^. For sufficiently large values of \n\, a reads 

m 

e^ , 

CJ(«) = -Lv^- (9) 

Going back to equation Q, we notice that the effective Hamiltonian of pristine graphene 
is invariant under a large number of symmetry operations in the isospin spaces besides the 
usual spatial translation and rotation invariances fTT|. First, the Hamiltonian commutes with 
the operators 

Av = a^Tx, Ay = (7,Tv, A, = (7oTz, (10) 

which, together with Aq = (Jo'^o, form a group (here, the subscript denotes the identity 
operator). Second, the Hamiltonian is also time-reversal and chiral symmetric; the 
corresponding operators are denoted by Tq = OxTx and Co = O^To, respectively. Notice that 
time-reversal is also accompanied by an inversion in the sign of the linear momentum as well 
as a transposition of the Hamiltonian matrix (it basically takes K to K' and vice-verse). A 
total of fifteen distinct symmetry operations can be constructed by combining Aq .v y.^, 7b, and 
Cq. For instance, {ToAxy^ H (TqAx) =H and [CqAx)^^ H [CqAx) ~ H (the latter exchanges 
n band branches but preserves the energy spectrum). 

Disorder will not only break translation and rotation symmetries but also affect these 
other invariances. Thus, to describe static disorder in general, we add to the Hamiltonian in 
(|5]l a term of the form 

Hdis = l^VijaiTj. (11) 

ij 

Depending on which matrix elements {V,;} are nonzero, certain pseudospin symmetries are 
broken while others are preserved. For instance, a generic nonmagnetic disorder may break 
all symmetries but time-reversal. In fact, the study of the effect of disorder on the transport 
properties of two-dimensional massless fermions described by Hamiltonians of the form (|5]l 
actually precedes graphene, and some important facts have been known for quite some time 
within the contexts of the integer Quantum Hall effect and d-wave superconductivity. Let us 
discuss a few important cases. 

Chiral symmetry is associated to the block off-diagonal form of H in equation (|5]l. 
Therefore, any disorder that introduces diagonal matrix elements or even a finite chemical 
potential (moving the Fermi energy away from the Dirac point) will break this symmetry and, 
in most cases, strongly affect the electrical conductivity. However, whenever Voj = V^j = 0, 
quantum corrections to the conductivity cancel each other to all orders in the disorder strength 
and there is no localization at the Dirac point llT4lll5l[T6l . The conductivity in this case takes a 
finite value which may depend on the disorder strength and whether the remaining symmetries 
of the Hamiltonian are broken or not [Tt] T3l. Bond disorder due to lattice distortions such 
as ripples, random magnetic fields, and dislocations and other lattice defects represented by 
non-Abelian gauge fields [18] fall into such a class of disorder. 

When chiral symmetry is not present, two important situations arise. First, if the disorder 
is sufficiently smooth, varying significantly only over length scales larger than the lattice 
constant, there is very little mixing between the two valleys; the K and K' points remain 
essentially decoupled and V,^ = V,-, = 0. For some time in the literature, the fate of massless 
Dirac fermions in the presence of long-range scalar disorder was unclear (T9\ |20l . In fact, 
some issues appeared even in the clean limit. For instance, a calculation of the Kubo 
conductivity of non-disordered Dirac fermions in the so-called ;r-flux phase model disagreed 
with the Einstein conductivity calculated in the thermodynamic limit I.20J . Only recently an 
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explanation for this discrepancy was provided [21 1, one reason being that since the neutral 
point is actually a critical point, the value of the clean conductivity depends on the order of 
limits. This indicates that one should also be careful when considering claims of universality 
in the disordered case. However, as we will discuss below, it became widely accepted recently 
that localization due to quantum interference is also absent in this case. Second, when disorder 
has a strong short-range component and intervalley mixing is significant, the opposite occurs 
and quantum corrections tend to normalize the conductivity to zero, turning graphene into 
an Anderson insulator II22I l23l l24l l25l . In this case, as far as transport is concerned, charge 
carriers in graphene behave much like those in an ordinary two-dimensional electron gas. 

Finally, it is worth remarking that far away from the Dirac point, graphene is no longer 
described by the Hamiltonian (|5]l since the dispersion relation ceases to be conic. For instance, 
trigonal warping changes the role of quantum corrections to the conductivity and enhances 
localization 1,26 J . We return to this issue in section |5] when we address mesoscopic effects in 
graphene. 

3. Nature of disorder in graphene 

It is commonly accepted that graphene has very few lattice defects. Therefore, intrinsic 
disorder tends to be weak even in exfoliated samples. However, extrinsic disorder is invariably 
present and is basically dictated by the synthesis method and by way the graphene sheet is 
supported. Chemical contamination in the form of adsorbates can occur during lithographic 
processing of field-effect devices. The focused electron beam used in the lithographic process, 
and imaging with SEM also introduce disorder, although the its exact effect on the graphene 
lattice is still under debate f28,"271. In addition, when graphene is laid on a substrate, lattice 
distortions (the so-called ripples) can appear due to a tendency in graphene to conform to 
the roughness in the substrate surface [29 \. It has been recently suggested that wrinkles and 
other singularities can also be formed, causing lattice distortions f3Q \ . Finally, most insulating 
substrates used in graphene devices are oxides which are rather prone to charge traps. These 
traps can be located either in the bulk of the substrate (but not too far away) or at the interface 
between the substrate and the graphene sheet and are a common source of disorder 

For clean, suspended graphene systems, extrinsic disorder can be substantially reduced 
II311 1321 . Nevertheless, some amount of disorder is still present, as shear and strain created 
by contacts and scaffolds typically induce corrugation in the graphene sheet f33] which can 
create electron scattering [34, 35, 36, 37, 38, 39, 40|. Ripples induced by thermal motion 
are rather unlikely at low temperatures, but can be an important source of scattering at high 
temperatures in the form of flexural phonon modes. 

It was recognized early that screening in graphene is rather poor due to the low density 
of states near the neutrality point. Thus, at low doping, charge traps in the substrate as 
well as Coulomb impurities (e.g. charged adsorbates) act as long-range electron scatterers 
Bn l42l |43] l44l |451 . It was shown by several authors that Coulomb scattering (even when 
screened) leads to a linear increase of the conductivity with electron density, matching quite 
well the majority of the experimental data for field-effect graphene-based devices. However, 
there are also a few other scattering mechanisms that can also produce a linear dependence of 
the conductivity on carrier density. 

There is one additional type of disorder in graphene that, albeit static, does not quite 
follow equation (fTTT l. It is induced by neutral adsorbates, such as atomic hydrogen, which 
tends to bind covalently to carbon atoms in the graphene sheet and locally distort both the 
lattice and electronic structure ll46l l47l . The result is a mid-gap or resonant state near the Dirac 
point that acts as a short-range scatterer, similarly to (but not exactly like) vacancies. This type 
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of disorder also induces a dominant linear dependence of the conductivity on carrier density, 
modulated by weak logarithmic term P8ll49l[50ll511 . Scattering by ripples in the graphene 
sheet can also, in principle, lead to a linear carrier density dependence |52|, although this 
situation is considered unlikely since it requires free-standing, equilibrium fluctuations in the 
graphene sheet which are unlikely to dominate when a substrate is present. Finally, a finite 
amount of wrinkles may also induce the same linear dependence lf30l . 

There is some recent experimental evidence that resonant scatterers may play a 
significant role in limiting mobility in graphene sheets deposited on oxide substrates by 
mechanical exfoliation f5J','54|, but their microscopic nature is still unclear. In fact, at present, 
the dominant mechanism of electron scattering in graphene is still under debate for both 
suspended and non-suspended graphene. In particular, for the latter, despite some compelling 
evidence |55 1, the widespread view that Coulomb impurities are the most important scattering 
mechanism limiting mobility has been recently challenged [56|. 

Finally, we mention that disorder also plays a role in other characteristics of graphene 
devices, such as 1 //noise ll57l l58l l59l . 

4. Neutral graphene 

Graphene at the neutrality point has some remarkable properties. As mentioned earlier, the 
conductivity is finite even though, in the clean limit, the density of states vanishes. Early 
transport experiments found a conductivity minimum of the order of e^/li, which raised 
the possibility that this may actually be a universal feature of graphene. It took some time 
for this issue to be settled since several theories also pointed to some universality in the 
conductivity of Dirac fermions at the neutrality point. For instance, in the presence of 
certain types of chiral disorder such as random vector potentials, it was well known that 
the conductivity is essentially unaffected by disorder (as long as the disorder is not too 
strong), taking the universal value of 4-e^/nh 1201 160^. For other types of chiral disorder, the 
conductivity at the neutrality point is also of order e^/ /z 1 17 1 . In fact, more than two decades 
ago, Fradkin used a coherent potential approximation to show that two-dimensional zero- 
gap semiconductors have a universal Drude conductivity of order e^/h at zero energy when 
intervalley scattering is neglected 1 19| . Some years later, Lee used a self-consistent Born 
approximation (SCBA) to show that a similar result holds for fermions at the nodal points 
of dirty d-wave superconductors fET]. In an early theoretical work, Shon and Ando used the 
SCBA to show that for both short- and long-range impurity scattering, the conductivity in 
graphene is equal to 2e^/nh per spin degree of freedom ll62l (see also ifTTI ). Coincidentally, 
this is the same value of ballistic graphene when W/L ^ 1, as seen in section|2] 

As more experiments began to show that the conductivity minimum can in fact vary 
widely from sample to sample M63I . thus contradicting the idea of universality, it became clear 
on the theory side that the SCBA is inadequate to describe Dirac fermions near the neutrality 
point. First, it does not correctly incorporate quantum effects such as coherent scattering 
and localization. Second, it can hardly be justified since at the Dirac point there is no clear 
expansion parameter controlling the approximation (the carrier wavelength diverges at the 
neutrality point for a uniform system). Third, and perhaps more importantly, real systems 
are far from being macroscopically homogeneous and uniform ||64] |65] l66l l67l . Strong 
fluctuations in the background potential can occur at length scales comparable to the mean free 
path /. This is not easily taken into account in the SCBA or in any semiclassical formulation, 
although some recent attempts have tried to circumvent this problem ll68l[69l . 

The sample-dependent value of the conductivity minimum shows that the nature and the 
strength of the disorder plays a crucial role in determining the transport properties of graphene 
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and a careful theoretical analysis must be employed. 

Let us first consider non-chiral, short-ranged potential disorder. Scattering induced by 
localized defects, neutral impurities or adsorbates can transfer enough momentum to carrier 
so that coupling between valleys across the Brillouin zone (i.e. the K and K' points) occurs. 
Several authors have considered this regime theoretically within the renormalization group 
approach ll22l l23l l24l [TT^ . They have shown that, as in the case of standard two-dimensional 
metals with no spin-orbit coupling, coherent backscattering in graphene also leads to carrier 
localization when the samples are sufficiently large, namely, when the localization length 
/loc is smaller than the sample linear size L. Graphene in this case belongs to an orthogonal 
symmetry class |fT3] l22l |23] l24l |26l ; in the one-parameter scaling language, it has a beta 
function j3((7) = dlncr/dlnL that is always negative, as shown in figure |2ja). This result 
is supported by several numerical simulations 170] |7T] |72l. However, it appears that, 
experimentally, L is either too small or the dephasing length is shorter than Zioc (we note 
that, in two dimensions, the localization length can be exponentially larger than the classical 
mean free path). Thus, the most common effect of short-range disorder at the neutrality 
point is a reduction (but not a complete suppression) of the conductivity. However, when 
short-range disorder is strong, variable-range hoping conductivity is observed, indicating 
localization. We will get back to short-range disorder when we consider mesoscopic effects 
and doped graphene. 
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Figure 2: One-parameter conductivity scaling functions proposed for graphene in the presence 
of (a) short-range and (b) long-range scalar disorder In (b): the dashed curve represents the 
standard symplectic class with an unstable fixed point at as, the dotted curve was conjectured 
in ||75l and the soUd curve corresponds to the result of II76II79I which is backed by numerical 
simulations II78II80|[8T1[821. The arrows indicate the direction of the renormalization flow. 



When the potential disorder is smooth at the atomic scale, intervalley scattering (and 
therefore backscattering) is suppressed. As a result, carriers can be described by a single- 
valley Dirac Hamiltonian with a random scalar potential ||20j|73J, 

H = vpip-a) + V{r)Go. (12) 

This Hamiltonian is symmetric under a pseudo time-reversal symmetry: H = GyH^Gy (notice 
that the real time-reversal symmetry operation would connect K and K' points instead). This 
Hamiltonian belongs to a symplectic class ifTS] l20l l22l |23] l24l l26l . It is well known that 
normal metals with spin-orbit coupling fall in this class and show a metal-insulator transition 
II74I . therefore, it seems plausible that the same could occur for graphene. In equation ( 112b . 
the pseudopsin (due to the sublattice structure) plays a role similar to the real spin in normal 
metal and is coupled to the orbital motion through the p ■ a term. From a theoretical viewpoint, 
until recently it was not clear whether a transition also exists for graphene in this class. Since 
it is difficult to fabricate samples with exclusively long-range disorder, experiments did not 
provide any insight into this question. 



Disorder and electronic transport in graphene 



8 



A field-theoretical treatment of this problem was proposed by two groups 1113117511761 . 
In particular, Ostrovky et al have argued that the appropriate non-linear sigma model includes 
a topological term which causes the appearance of a novel stable fixed point at a finite 
conductivity a* in addition to the usual symplectic unstable fixed point Os w 1 Ae^/h 17511771 . 
The suppression of localization was also obtained by Ryu et al ll76l using a combination of 
analytical and numerical calculations. These authors also derived a topological term which is 
explicitly constrained to be either or ;r by symmetry. In contrast to |75|, their treatment of 
the global anomaly encoded in the topological term yields only one stable metallic fixed point 
and no transition. Both results are illustrated in figure|2|b) using a one-parameter scaling beta 
function. 

The derivation of the topological term in IfTSl l75l indicating the validity of the single- 
parameter scaling of the conductivity, the existence of a stable fixed point, and the consequent 
metal-insulator transition requires a non-perturbative investigation. This motivated several 
numerical studies that employed a variety of techniques ifTS] 1791 [80l 1811 l82l . Only a strictly 
positive beta function was obtained and no evidence of a new fixed point was found, thus 
corroborating the prediction of Ref. |76|. This supports the idea that the main role of the 
topological term is just to suppress the conventional symplectic metalinsulator transition. 
In all studies, the average conductivity followed a simple scaling law, increasing with the 
logarithm of the system size such that j3 (a) = a/cT. While one could in principle expect the 
constant a to be universal, the situation is somewhat more complex, particularly for finite-size 
systems (with the exception of |79 1), since simulations are carried out at zero temperature and 
therefore in the fully coherent limit. In addition, most simulations obtain the conductivity 
by computing first the two-terminal conductance G and can be influenced by contacts and 
geometry. 
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Figure 3: Average conductivity of a graphene strip (armchair edges) as a function of rescaled 
length for three different aspect ratios. Several values of disorder strength {Kq = 1,2,4,8, 16) 
and correlation length (^ /a — 2,4,8) were used. The number of realizations range from 200 
to 1000. The dashed line is a visual fitting to the large-L portion of the numerical data. The 
dotted lines indicate the ballistic regime, equation (|7]|. 
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In our numerical simulations, we have found that the constant a deviates from the 
universal value of e^/nh expected for the symplectic class. We use a single-band tight- 
binding model on a honeycomb lattice with an added long-range Gaussian correlated potential 
disorder PH . namely, 

(VWV(/))^/.of^exp(-t_^), (13, 

where Kq parametrizes the disorder strength and is the disorder correlation length. Figure 
Oshows the finite-size scaling of the average conductivity using the recursive Green function 
method lISTl . In essence, this method mimics a two-terminal transport measurement. Three 
different aspect-ratios were considered as well as several values of and ^/a. The average 
conductivity obtained from the linear conductance is plotted versus L/l, the ratio between 
system length and the disorder mean free path /, which depends on Kq and t, . Two clear 
regimes can be identified. For L/l ^ 1, the probability of an electron being scattered by 
disorder as it traverses the sample is very small. This corresponds to the ballistic regime, 
where scattering occurs mainly at the sample edges and transport properties are dominated 
by the sample geometry. In contrast, when Ljl^ 1, the system becomes diffusive. Now, 
electrons are scattered multiple times by disorder in their path across the sample. Scattering 
at the sample edges is no longer dominant and transport depends weakly on the geometry. 
Figure 13 clearly shows this crossover In the L/l Q limit, the simulations approach the 
analytical result (indicated by arrows) for clean graphene as given by equation For the 
diffusive regime, L/l^ 1, the conductivity is proportional to ln(L//), in agreement with the 
non-linear sigma model prediction II17I . The mismatch between the numerical prefactor for 
the logarithm and the value characteristic of the symplectic class may be related to the finite 
contact resistance. 

These simulations suggest an explanation for the transport experiments at the charge 
neutrality point. In the coherent diffusive regime, the conductivity has significant sample-to- 
sample fluctuations and its average shows a weak (logarithmic) dependence on the mean free 
path. Typical diffusive experimental samples have L// w 1 — 10 and a ~ 4e^//i, similarly to 
what is seen in figure [3j 

An intuitive, non-rigorous, understanding of the lack of localization in the presence 
of long-range disorder is possible by invoking the so-called Klein tunneling f83] [SH [85l . 
In the absence of intervalley scattering, as carriers move across charge puddles induced by 
long-range potential fluctuations, chirality prevents backscattering. When the Fermi energy 
cuts across a potential barrier, carriers are converted from particle to hole (or vice-verse). 
Chirality conservation then requires a forward moving particle hitting the potential barrier at 
a normal angle to be scattered as backward moving hole (interband tunneling). In this case 
there is perfect charge transmission, as the real electron involved in the process continues 
to move in the same direction. For other incidence angles, a finite probability of particle 
backscattering appears and it is actually a periodic function of the angle. This effect has 
been recently observed experimentally |'86l. Klein tunneling through charge puddles (p- and 
n-doped regions) is the basis of an attempt to explain the conductivity in graphene near the 
neutrality point using a random resistor network model lISTl . 

5. Mesoscopic effects 

Mesoscopic corrections to the conductivity have been studied for more than 20 years in 
numerous systems, such as metals and two-dimensional semiconductor heterostructures 
at low temperatures ll88l [89l l90l 1911 l92l . Weak localization and universal conductance 
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fluctuations are the most ubiquitous manifestations of mesoscopic effects in transport. The 
weak localization correction to the conductivity can be understood semiclassically in terms 
of quantum interference between self-crossing paths originated by disorder scattering. In the 
loops, an electron can propagate in the clockwise and counter-clockwise directions and these 
two trajectories interfere. The effect is maximal when the system is time-reversal symmetric; 
for instance, in normal metals with negligible spin-orbit coupling, it leads to constructive 
interference and a decrease in the probability of carriers to move forward. This interference 
effect gradually disappears as time-reversal symmetry is broken, for instance, by an external 
magnetic field. In diffusive systems, mesoscopic effects are also responsible for sample- 
to-sample fluctuations in the conductance, 5G, of order e^/h, irrespective of the magnitude 
of G itself I193I . This phenomenon is known as universal conductance fluctuations (UCF). 
The magnitude of quantum interference effects strongly depends on £^ //, the ratio between 
the electron dephasing length and a sample specific length scale, such as the electron elastic 
mean free path or the linear size of the sample L. When iij, ^ L, this magnitude saturates 
to a constant value; in convention metals, this value is solely determined by the symmetry 
class and the sample geometry. Since, for any given sample, £^ decreases very rapidly with 
increasing temperature, low temperatures are usually required to observe mesoscopic effects. 

The effects of quantum interference in the electronic transport in a honeycomb lattice had 
been theoretically addressed |94| even before graphene was first synthesized. Suzuura and 
Ando focused their analysis on the weak localization peak, namely, the enhanced resistivity at 
zero magnetic field. Using the standard disorder diagrammatic approach for diffusive systems, 
which amounts to a perturbative expansion in powers of (fcp/)^', it was shown that the weak 
localization correction depends strongly on the spatial range of the disorder (Notice that this 
expansion is only formally justifiable away from the charge neutrality point, when kp is finite.) 
As mentioned in section |2l due to chirality, backward scattering in graphene is suppressed 
for long-range disorder. In this case, electrons in the vicinity of distinct valley points do 
not couple and the system belongs to symplectic class, leading to an anti-weak localization 
correction (negative magnetoconductance), namely, a valley rather than a peak in the curve 
of resistivity versus magnetic field. Short-range (non-magnetic) disorder, in contrast, can 
break all symmetries except the time-reversal one. In this situation, the system belongs to the 
orthogonal class and there is weak localization (positive magnetoconductance). 

The theory was later extended to include trigonal warping effects li26l . The diagrammatic 
expression for the weak localization correction, Aa{B) = g{B) — cj(0), as a function of 
magnetic field reads 



-2% 



(14) 



where F{z) = lnz + 'P(l/2+ 1/z) and ^ is the digamma function. The magnetic field is 
cast in terms of ' = AeDB/Ti, where D = v'^Tti/2 denotes the diffusion constant. The three 
other times scales correspond to: the dephasing time T^, which is related to the dephasing 
length as t.^ = (Dt^)'/^; the intervalley scattering time Ti, which is due to sharp, atomic-like 
disorder features; and the intra-valley scattering time T*, which is related to extended defects, 
dislocations, ripples, etc.. The derivation of this expression assumes that the momentum 
relaxation time, Zp, is the shortest of all time scales and is caused by charge impurities (and 
not by atomically sharp defects), so that it does not affect chirality. 

The first experimental report on mesoscopic interference effects in exfoliated graphene 
observed a strong suppression of the weak localization peak, as compared to the typical value 
of e^ /h, independently of doping even at very low temperatures |95 |. This behavior was 
attributed to the presence of ripples in the graphene sheet: the effect of corrugations in the 
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electronic structure can be translated into an effective random gauge field which destroys 
quantum interference in the same way as a random magnetic field. Further experiments in 
different exfoliated graphene samples f96l, where temperature and carrier concentration were 
systematically varied, observed much larger weak localization peaks. Experiments in epitaxial 
graphene, on the other hand, showed an extremely sharp anti-localization peak 11971 . 

These observations are consistent with the theory of |26 | when the time scales appearing 
in equation (114b are properly adjusted. In line with the qualitative discussion presented above, 
a negative magnetoconductance correction or weak anti-localization. Act < 0, corresponds to 
the case when T, — > °°. Weak localization and positive magnetoconductance. Act > 0, occurs 
when T,- 0. 

More interestingly, the theory predicts a region in the parameter space where a transition 
from positive to negative magnetoconductance corrections occurs. To see that, it is convenient 
to consider a weak magnetic field, such that one can approximate the F{z) function by its 
leading order expansion, F[z-^0)~ Then, 

, / _ „T \ 2 . 

1 2 



' 3 /i I Oo , 



1 



(15) 



where (S>Q — h/e is the flux quantum. Based on this expression, Tikhonenko et al |98 1 proposed 
a diagram reproduced in figure |4] By following the dashed blue line, one can cross from a 
localization to an anti-localization region. Since the scattering times T; and are essentially 
temperature independent but is not, one can use temperature to move along this line. 
Another way to move in this parameter space is to vary the carrier density: low density and 
high resistivity is associated to short scattering times, while high densities and low resistivity 
goes in the opposite direction. The experimental data 1981 nicely support this analysis. 





localization 


anti-localization ' 





Figure 4: Diagram illustrating the weak localization-weak anti-localization crossover 
Adapted from |98|. 



The picture that emerges from this analysis is that transport in exfoliated graphene on 
Si02 substrates is influenced by both long- and short-range scattering processes. For standard 
experimental samples, the strength of short-range scattering processes is such that hoc ^ L. 
Recent experiments ['99' '1001 show that the conductivity is strongly suppressed when the 
concentration of atomistically sharp defects is artificially increased, likely because l\oc ^ L. 

A quantitative study of mesoscopic quantum corrections to the conductivity at the charge 
neutrality (Dirac) point lies outside the validity range of the diagrammatic approach. The 
interplay between localization length and the finite system length indicates that numerical 
tools are necessary to gain further theoretical insight in this case. Qualitatively, however, the 
main features are similar to those obtained at finite doping since, particularly when charge 
puddles are present, since in this case A:p is nonzero locally. In figure |5] we show the result 
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of numerical calculations of the magnetoconductance for a strip with armchair edges, at the 
neutrality point and away, for both short- and long-range disorder. The calculations were 
performed with a tight-binding model and a Gaussian disordered potential, similarly to those 
presented in figure |3] Notice the transition from weak localization to anti-localization as the 
range of the potential is made larger than the lattice spacing. 
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Figure 5: Average magnetoconductance of an armchair strip (W/L = 3) at (Vgate = 0) and 
away from (Vgate > 0) the neutrality point for short- — Q.5a) and long-range (£, — 4.0a) 
disordered Gaussian potentials. <1> is the total magnetic flux perpendicular to the graphene 
sheet, <Po = h/c is the flux quantum, and f is the tight-binding nearest-neighbor hoping 
amplitude. A total of 500 reahzations for each case have been used in computing the averages. 



Quantum interference among the multiple paths that carriers can take as they move across 
the sample also cause conductance fluctuations. Despite their random, aperiodic nature, these 
fluctuations are reproducible when an external parameter is swept back and forth. They are in 
essence quantum fingerprints of the microscopic configuration of the disorder in the sample 
and are usually quite evident when L ^ A proper way to quantify UCF is to evaluate 
the autocorrelation function C{AX) — {5G{X + AX) dG{X)), where X represents an external 
parameter such as a magnetic field and 5G is the deviation of the conductance from its average 
value. 

Early transport measurements in graphene were able to identify mesoscopic fluctuations 
moil [95] 11021 11031 . but only recently more careful quantitative analyzes were carried out 
by several groups [104, 105, 106, 107] 11081 . On the theory side, numerical simulations 
showed that long-range disorder in graphene can produce anomalously large conductance 
fluctuations flOQ], while for short-range disorder the magnitude of the fluctuations is similar 
to that observed for normal diffusive metals. This result was explained within the framework 
of the diagrammatic perturbation theory by Kharitonov and Efetov 01 101 . who highlighted 
the crucial role played by valley symmetry breaking in determining the magnitude of the 
conductance variance varG — C(0). For Dirac fermions, three possibilities arise: when all 
contributions to scattering are relatively weak, varG is four times larger than the standard 
value for normal diffusive metals. When intravalley scattering is strong or triagonal warping 
is present (high doping), varG is enhanced by just a factor of two. Finally, when intervalley 
scattering is strong, varG takes the standard value. They also showed that when time-reversal 
symmetry is broken by an external magnetic field, varG is further reduced by a factor of 
two, as it is the case for conventional metals. Similar results were also obtained in 1111], 
where the role played by geometry (quasi-one-dimensional versus two-dimensional) was also 
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emphasized. However, as was the case for the diagrammatic calculations of weak localization, 
both theories are only applicable away from the neutrality point, as they rely on kpl ^ 1 . 

Experimentally, there is still some disagreement about the magnitude and other 
characteristics of UCF in graphene. One group |104| reported that UCF are suppressed near 
the neutrality point in both bilayer and trilayer graphene. The same effect was observed in 
monolayer samples by the Illinois group |108|, who related it to a marked reduction of the 
coherence length £(j, near that point. However, another experiment on monolayer and bilayer 
samples saw the opposite behavior II 1071 when mesoscopic fluctuations were collected at a 
fixed value of an external magnetic field, while fluctuations collected by varying the magnetic 
field were roughly independent of the gate voltage. The authors in |107| argued that this 
discrepancy is due to the dominant presence of strong charge inhomogeneities (puddles) 
in their samples: while varying an external magnetic fields affect little the distribution of 
the puddles, gate voltage variations are much more effective in scrambling the conductance. 
There is also some discrepancy among the effect of magnetic fields on the UCF amplitude: 
while 1105 1 sees the expected factor of two, in 01 071 this only occurs near the neutrality point, 
while in fl08 | no reduction is observed. 

The temperature dependence of the UCF amplitude was measured by several groups. 
In particular, |105| finds a clear inverse power-law behavior at high doping, typical of 
conventional diffusive metals, combined with a saturation at low temperatures. However, the 
authors in |106 | observed an exponential suppression of UCF with increasing temperature. 
While the former is backed by the standard theory of dephasing in metallic systems LI 12j . the 
latter has yet to be understood. 

Finally, we mention a recent experiment where both weak localization and UCF were 
observed when an in-plane magnetic field was applied to graphene on Si02 1113]. The 
appearance of quantum interference in this case can be directly connected to rippling in the 
graphene sheet. 

6. Doped graphene 

Doped graphene shows remarkable high field-effect mobilities, even at room temperatures 
111 14, ,1 151 . This is a key property that makes this material a potential candidate for future 
carbon-based electronic devices. The dominant source of scattering limiting the mobility jj, in 
graphene is still under debate. As already mentioned, at low temperatures, it is believed 
that charge impurities [41. jl 161 |69l Ill7|, substrate effects, resonant scattering |48, 49J, 
corrugations 11521 . and strain are the main sources of disorder. 

With increasing doping, the conductivity of graphene depends only weakly on 
temperature and grows almost linearly with carrier density, in contrast with a \n\^^^ 
dependence predicted for clean graphene, equation (|9]l. These features were first explained 
by considering charge impurity disorder, likely due to charge trapped in the substrate 
BTI II 161 l69l II 17L These theories also predict that the mobility should be inversely 
proportional to the impurity concentration nimp (for uncorrected scatterers) and to strongly 
depend on the dielectric constant. So far, there is no experimental consensus on these 
properties. While experiments using potassium ions provide evidence that jj, oc I /n^^p \ 118], 
experiments with gaseous adsorbates have only showed a weak dependence of jJ. with «in,p 
111 191 . There is also an experimental controversy regarding the effect of the dielectric 
constant on the conductivity: The Manchester group measured only a modest change in 
by immersing graphene in high-fc environments, such as ethanol and water ll56l . while the 
Maryland group observed a significant change in fi by adding just a few monolayers of 
ice to a graphene sheet 1,55 J . Recently, experiments with hydrogenated graphene revealed 



Disorder and electronic transport in graphene 



14 



a strong short-range scattering component to the mobility, apparently more dominant than 
charge impurities in some samples M53II54I . 

Let us now discuss the main theoretical approaches to address the mobility in graphene. 
For finite doping such that kpl^ 1, it is possible to calculate the conductivity within much a 
simpler theoretical framework than that used at zero doping, where one deals with a non- 
perturbative problem (see section |4|i. In addition, as one moves away from the charge 
neutrality point, mesoscopic corrections to the conductivity lose their importance, since 
5a/a < 1. 

Assuming the system to be homogeneous and diffusive, I, the simplest approach to 
describe the linear transport properties of graphene is given by the Boltzmann theory |l5]|6]|7l, 
which is widely used as a guide for interpreting experimental results in doped graphene. This 
semiclassical approach gives the conductivity in terms of the transport scattering time Tq-, 
namely 

a='^ j dE(^-^v{E)vlx^{E), (16) 

where v(£') is the density of states, f{E) is the Fermi-Dirac distribution function, and Tq- is 
usually calculated using Fermi's golden rule, namely, 

^ - ^/ (l-k-l^)F,^3{Eik) -E{k)), (17) 

where q= — = 2A:sin(0/2) is the momentum transferred, V{q) is the Fourier transform 
of the scattering impurity potential, and k k = cos Q. Chirality effects are captured by Fij^^, 
which is F^t = (1 + cos 0)/2 for graphene. 

By modeling short-range disorder with zero-range delta-like scatters V{r) = mq 5 (r — 
r,) ll62l . one readily obtains that T^' °^ ^(s)- As a result, the conductivity does not depend on 
the chemical potential E^ or on the carrier density n, which is obviously incompatible with the 
experimental results. This observation is frequently used to dismiss the importance of short- 
range scattering in graphene. Thus, one has to bear in mind the limited applicability of this 
model. At high doping and considering any kind of realistic of disorder with a spatial range 
(flange, such that ^p^range 7^ 0, modeling disorder by zero-range scatterers is hardly justifiable. 
This model is also not very helpful for small doping either, since the Boltzmann semiclassical 
approach breaks down for small k^l. Even more importantly, it was realized very early P1201 
that short- (and nonzero-) range disorder leads to resonant scattering. In this situation, the 
Fermi golden rule, equation ( 117b . breaks down and Ttr should be calculated instead as follows 
01211 : (i) use resonant scattering to calculate the transport cross section Oj, as it has been 
done, for instance, in [6 | for a model where randomly placed hard discs of radius R are the 
disorder source, (ii) use 1 /t^ = n[m-pV-pOT{kR) II121I instead of equation ( fTTI ). By proceeding 
this way, it was shown [i^lllSlllgl that resonant scatterers give raise to mid-gap states and lead 
to a conductivity quasi-linear in n, depending on the radius R of the scattering discs. 

The Boltzmann transport theory is customarily used to address the effects of charge 
impurities in graphene 111 161 1411 l44l . In this case, the disorder scattering potential V{q) 
becomes Vc{q) / £{q), where Vc{q) = Ine^jiKq) is the Fourier transform of the Coulomb 
potential in an effective background dielectric constant fc, whereas e(g') is the graphene 
dielectric function. Within this disorder model, the influence of the substrate on the electronic 
transport can be directly accessed by properly modifying fc. This is the theoretical basis for 
the experimental work in lf55l l56l . The dielectric function 

e{qJ)^\^Vc[q)X\{qJ) (18) 
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with Yl{q,T) denoting the irreducible polarizability function given by the standard bare bubble 
diagram II411 11221 l44l . The non-analytical behavior of n(^) at ^ 2k f leads to Friedel-hke 
oscillations in the screened potential in real space, namely, 0(r) ^ co^{2kpr)/r^ ||84l 11221 . 
Both n(^,r) and the chemical potential pi{T) temperature dependencies have a strong 
influence on the conductivity. For low temperatures, such that T/Jp ^ 1, the conductivity 
depends very weakly on T, whereas for T/Tp ^ 1 one finds [123] that o °^ {T /Tp)^. 

The Boltzmann theory has been also used to calculate the conductivity of graphene in the 
presence of static ripples [49 1 and phonons P8lll24lll25l . with limited success in explaining 
the experiments, in particular the temperature dependence of a. Let us also mention that 
the standard Boltzmann approach was used to calculate the graphene conductivity for the 
correlated Gaussian disorder model discussed in sectionlH yielding 11261 



2 r 



= ^ {2n\n\^y/' + 0{\n\er- 



(19) 



Since Kq oc 111091 . the conductivity increases linearly with the inverse disorder 

concentration n^^. For high caiTier concentrations, such that nt^^ 3> 1, C7 is proportional 
to 

The next standard approximation level beyond the Boltzmann transport theory is the self- 
consistent Born approximation, as mentioned in section |4] The starting point for the latter is 
the Kubo formula, which gives the conductivity in terms of Green's functions. While the 
standard Boltzmann approximation treats the scattering within the Fermi Golden rule, the 
SCBA prescribes an efficient way to encode the main scattering processes of an electron in a 
disordered potential into the Green's function self-energy (see for instance 01271 ). 

Ostrovsky et al (Tl] presented a thorough study of the conductivity due to generic 
(long- and short-range) Gaussian disorder in monolayer graphene at finite bias using the self- 
consistent Born approximation and a renormalization group (RG) analysis ll23l . While for 
non-resonant (or weak) impurities the SCBA and the RG analysis give similar results, for 
resonant scattering the RG conductivity renders 

. ^ ^2 \n\ 2 A2 

where A is a momentum cutoff. This result is very similar to the one obtained for the hard disk 
model |6 1. Note that the RG analysis for Gaussian disorder, including resonance scattering, is 
significantly different than the one obtained using Fermi's golden rule, equation il9i 1 126 1. 

The RG analysis i23l[T7i has been also recently employed to study the conductivity in the 
presence of randomly placed adsorbates 1 128 1. For resonant states far from the Fermi energy, 
it predicts quite distinct and asymmetric curves for a{n). For instance, a(n) can shown an 
almost insulating-like behavior for p-type graphene and be metallic-like for n-type graphene. 
This is consistent with some, but not all, experiments with adsorbates. 

Before concluding this section, it is worth mentioning that although typical, good quality, 
micrometer-size graphene samples show L/l > 1, many do not belong to the regime L/l ^ I. 
Hence, deviations from the diffusive theory presented above, such as in [63], should be 
expected. Numerical simulations provide insight into the ballistic-diffusive crossover regime 
at finite doping 181]: with increasing disorder strength, the conductivity dependence on the 
carrier concentration moves from a sublinear dependence, resembling that of equation (|9]l, to 
a superlinear dependence. 
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7. Conclusions and outlook 

In this review we attempted to highlight the most important developments related to the 
role of disorder in the electronic transport properties of graphene. We emphasized the role 
played by symmetry breaking and discussed the various types of disorder commonly found 
in exfoliated graphene sheets. Although some aspects of the interplay between disorder and 
quantum coherence still require theoretical and experimental investigation, much about strong 
and weak localization and mesoscopic fluctuations of conductance in graphene is already well 
understood. However, there is still no consensus about which scattering mechanism plays 
the dominant role in limiting mobility in both suspended and non-suspended graphene field- 
effect devices. The widespread view that charge traps in the substrate are the most effective 
scatterers of carriers in high-mobility graphene on oxide substrates has been challenged by 
recent experiments. 

Several important topics were not covered in this review. For instance, disorder 
(intrinsic or extrinsic) plays an important role in inducing spin-orbit coupling in graphene 
and, consequently, in reducing spin relaxation times ||129II130II131|[T321 . Thus, investigating 
how scattering mechanisms of various forms affect spin relaxation is critical for advancing 
graphene as a spintronics material. Another situation where disorder is critical is in 
nanoribbons 013311134)11351 . In general, it is fair to say that there are still many open problems 
to explore in electronic transport in disordered graphene systems. 
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